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Preface 


The work tle,scril)cd in this report was performed hy the Teleeomnuinications 
Division of the Jet Propulsion l^dmralory. 
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Abstract 

Radio l)cacons with distinguishing signatures exist in nature as pulsating radio 
sources (pulsars). These objects radiate well determined pulse traii-s over hundreds 
of megahertz of bandwidth at radio frequencies. Since they are at known posit ions* 
they can also l>e used as navigation beacons in interplanetary space. Pulsar signals 
are weak and dispersive when viewed from earth. If an omnidirectional antenna is 
connected to a wideband receiver (200 MHz bandwidth centered at 2(X) MHz) in 
which dispersion effects are removed, nominal spacecraft position errors of 1500 
km can l>e obtained after 24 h of signal integration. An antenna gain of 10 dB would 
produce errors as low as 150 km. Since the spacecraft position is determined from 
the measurement of the phase of a periodic signal, ambiguities occur in the 
position measurement. Simultaneous use of current spacecraft navigation .schemes 
eliminates these am[)iguities. 


JPl TECHNICAL REPORT 32-1594 


Interplanetary Navigation Using Pulsating 
Radio Sources 


I. Introduction 

A .scientific mission to the outer planets will result in an 
.'m proved epbemeris of these pUincts. Sul)sequent missions 
would then become easier, hence cheaper, l>ecause less 
maneuvering would be recjuired. Imagine a spacecraft 
located in a spherical coordinate system centered at the 
sun. The ranj^e of a particular spacecraft can l>e measured 
with uncertainties as low as 5 m. Angular positions of the 
spacecraft are relatively larger and increase with increas- 
ing range. For example, an error of 0.1 arc-s in angular 
position is an error of 3(X) km at the distance of Jupiter 
and 23(K) km at the distance of it'luto. These errors are 
small when compared to thr ephemeris uncertaiiities. In 
late 1973, these uncertainties were 8(X) km for Jupiter and 
^X),(XX) km for Pluto. It is then clear that a flyby of a 
spacecraft will lead to l>etter ephemerides. This report 
discusses the u.se of naturally occurring pulsating radio 
sources to fix a spacecrafts position within a given 
coordinate system. 

One is reminded of the Lor;;n-C navigation scheme, in 
v'hich the navigator can locate his position on the surface 
of the earth by obtaining relative time delays between a 
master radio transmitter and one or more of its slave 
transmitteis. He then locates himself at the intersection of 


two circles. A similar system could he used for interplane- 
taiy navigation by placing radio transmitters on or near 
selected planetary surfaces. 

A natural navigation system is available to us in the 
form of pulsating radio sources (pulsars). When pulsars 
were first discovered, the possibility that they were the 
beacons of an intragalactic navigation scheme was 
considered. (This was known as the LGM hypothesis, 
where LGM = little green men.) This hypothesis was 
soon abandoned when it was found that pulsars radiate 
over many hundreds of megahertz. In this discussion, we 
turn this hypothesis around and ask how we might use 
such objects. As in the case of the Loran-C transmitters, 
these sources are at known locations, and they emit pulses 
of radiation a few milliseconds wide at regular intervals 
(typically 0.5 s) and at predictable times. They are within 
our galaxy but outside the solar system. They are therefore 
l>eyond o”r control, but have useful lifetimes of several 
million years. They are intrinsically powerful emitters, 
averaging 10*^ J in each pulse between 100 and 10,000 
MHz (that is, the output of Hoover Dam for 3 years). 
However, the distances to these sources are so large (up to 
10^2 m) that energy densities in the wavefronts are often 
attenuated more than 450 dB. Radiometers on or near 
earth will then detect pulses with only 10-2^ J/m^ of 
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antenna area, su^'eslinj' that long integration times will 
l)c needed to obtain usable signal-to-noise ratios. 

The measurement problem is defined in Section II in 
only enough detail to allow a calculation of expected 
errors in the spacecraft position. In Section III, we 
consider the appropriate properties of pulsars and their 
surroundings. An analysis of a navigation system using 
pulsars as the l>eacons is pre.sented in Section IV, and 
concluding remarks are made in Section V. 

II. A Definition of the Navigation Problem 

The position of the spacecraft in a particular coordinate 
system is calculated by comparing the arrival time of a 
particular pulse at the spacecraft with the arrival time of 
that pulse at the origin. Alternatively, we mea.sure the 
pha.se of the pulse train at the spacecraft relative to a 
given epoch and compare it with a similar mea.surement 
at the origin of the coordinate system. The mea.surement 
process involves errors l>ecause of noise fluctuations from 
the receiver, the background radiation from the galaxy, 
and the pulsar itself. The size and orientation of the 
ellip.soid describing the uncertainty the spacecraft 
position can be estimated if enough system parameters are 
known. 

A, Measurement of Position 

Consider a pulsar located in the direction given by unit 
vector Up and a spacecraft whose position is defined by a 


Z (S = 90') 



Fig. 1 Barycentric coordinati, systam B containing the pulsar 
and the spacecraft Sp. The objective is to n>easure the 
components of the spacecraft vector S 


vector S (Fig. I). The difference in the tiinc-of-arrival 
(TO A) of a pulse or wavefront at the spacecraft and at 
origin B ' » 

7 - - S • Up (I) 

where c is the velocity of propagation, T is the (plant it v 
measured, and the components and of S arc the 

unknowns. A minimum of three independent eipiations are 
needed, so at least three pulsars must l)e included in E(j. 
(1). We therefore make the following measurements: 

7, == “ (S, cos (6.) cos (a,) 
r 

+ Sy cos (6.) sin (a,) + S, sin (6,)1 (2) 

where the subscript refers to a particular pulsar and takes 
on values of I. 2. and 3. The angular coordinates and 
»re the right ascension and declination, respectively, of 
the ith pulsar. Equation (2) can l)e written in matrix form 
as 

{n =-(/4XS) (3) 

C 

where (T) and (S) are 3X1 matrices and (A) is a 3 X 2 
matrix. We then solve for (S) if the inverse (A *) of (A) 
exists: 

/sA 

Sy = (S) * (^A‘*)(7> (4) 

VJ 

B. Measurement of Time-of-Arrival 

The errors assoc iated with the measurements descril)cd 
hy E(|. (4) depend on how the mca.su rements arc made. 
The mea.surement procedure is represented in Fig. 2. In 
this analysis, we a.s.sume that the signal, hand-limited with 
a rectangular pa.s.sband, is .s(juare-law delected so that the 
receiver output is proportional to the average power. The 
detected signal is then passed through a low-pass filter to 
reduce noi.se fluctuations. The filtered signal is filtered 
further by correlating the output with a noise-free 
template representing the average shape of the received 
pulse. That is, the signal is pas.scd through a matched filter 
to maximize the signal-to-noisc ratio (S/N). It is assumed 
that the average pulse shape or template has l)ccn 
obtained in a separate observing program. It is the result 
of adding several thousand pulses together and is 
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Fig. 2 R«pr«Mntation of the moaMiremtnt of tima*of- arrival. Tho trror in tha moasuronfont is duo 
to noiso proportional to (recoivor noise T|^ and sky noito Tp) within bandwidth AF 


represented by /(/) in Fig. 2. The correlation process is 
performed such that the time delay t l>etween the 
template and the signal is measured relative to a 
particular epoch /q- The maximum value of the correlation 
function c(t) occurs at delay The delay iu the TO A 
of the pulse or the phase of the pulse train relative to f<). 

Random fluctuations in the detected signal cause 
random fluctuations in (ir), so there is an uncertainty 
associated with the measurement of t„. The rectangular 
RF filter has a bandwidth of AF. Consider the integrator 
following the square-law detec tor to l)e a rectangular low- 
pass filter with a bandwidth of AF/ . The dc output of the 
detector is proportional to TpAF* where Tp is the antenna 
temperature of the pulsar. The fluctuating component of 
the detector output has an autocorrelation function R„(/). 
Represent *he pulse template by whe^e the maximum 
value is 1.0. It can then be shown that the variance of a 
TO A measurement is (Ref. 1) 

. . 4L ,5, 

'* “ SUHt 

where T, is the equivalent noise temperature of the 
system. We have defined an integration time At * 

1/2AF/ . (See Ref. 2 for a discussion of integration times.) 

Equation (5) is valid if is less than the pulse width. 

It seems rea.sonable that we can define time more 
accurately with a strong, narrow pulse than with a weak, 
broad pulse. The.se facts are demonstrated in the two 
factors of Eq. (5). 

C. Eriors in Position 

Three particular pulsars are included in the navigation 
scheme. If a large number of measurements were 
performed of each T| represented in Eq. (3), the errors 
would be independent and distributed normally about 


some tnie value <T,>. The covariance matrix of the series 
of mea.surements is then diagonal and written a' 


A 


CTi^ 0 0 

0 CT2^ 0 
0 0 a3^ 


( 6 ) 


where the subscript on has l>een simplified. The joint 
probability density function of the deviations from T, is 
then the normal distribution 


(7) 


where ^ <T„>. 


The errors in .snacecraft loc'ation can l>e found from the 
coordinate transformation of Ecj. (4). The spacecraft 
positions and will l>e normally distributed almut 

true mean values S^>, and The random deviations 
rj “ S, - etc., will have a covariance matrix (/i) given 
by 


(ft) * (A'‘)(A)(A-*)’‘ (8) 

(T denotes transposition of rows and columns). The 
probability density function of the errors is then 

/Kri. rj, rj) - [~ 

(9) 

where n, m ■ 1, 2, 3 and denotes a particular 

element of the symmetric inverse matrix (ft'O. 

The suir.ination in Eq. (9) is a <|uadric form F^ « 
(R)’(fi*M(R), where (R) is the column matrix of errors r^, 
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vhich \fi : ilffect a 





ro, and r^. The form is positive definite, so a partieular 
value of F^^ defines an ellipsoid of a certain size centered 
at the Icxatiun calculaled in Kq. [i). Also, the particular 
' ahie of determines the probability that the spacecraft 
is actually located within the ell*;>soid. For example, when 
= 4 the prohahility that the vector lies within the 
ellipsoid centered on S is 0.95. It is then desirable to view 
the error ellipsoid in a coordinate system aligned with its 
principal axes. This is equivalent to diagonalizing the 
matrix (/t *). If the new inverse covariance matrix is 
then we need a coordinate transformation (^/, usually a 
series of rotations, such that 

{F-'y = 

If the columns of {(J) are !uadc up of the components of 
the eigenvectors of i/c M, then (()) will have the desired 
properties (see, for example, Ref. 3). 

Let Xj, X2, and X3 represent the eigenvalues of {/i ‘) and 
< n2’ ^’n.3^ represent the eigenvector associated 
with X„. The components of e„ are projections of e„ onto 
the original reference, frame defined in Fig. 1. The 
eigenvalues are found from finding the three roots of 

1(^-1) X(/)l = 0 (10) 

where (/) is the identity matrix. The eigenvectors can be 
found by noting that, letting (p) = to simplify 

notation, 

Pi 2^32 F Pl3^*13 ~ ^ (Ifa) 

P2I<’2I (p22 ^2)^'22 P23^'23 ® (Hb) 

P31<'3I P32^‘32 ^P33 “ ^3)^'33 “ ^ (Hc) 

E(]uations (lla-c) can be viewed as the expansion of the 

determinant l(p) ~ X„(/)l, which must equal zero. 
Therefore, the vector components can be the 

cofactors of the elements in the nth row of the matrix ip) 
“ X„(/). Normalizing the ve( tor c„ such that 

3 


produces a unit vector aligned along the principal axis 
a.ssociated with the half-size X„. 


III. Properties of Pulsar Signals 

Several important phenomena and facts ohseived in or 
about elect romi’gnelic waves from pulsars are presented 


and disciis.sed here. Only the proper ' 
navigation attempt are discussed. 

A. Pulse Period and Pulsar Location 

The interval between the pulses (pulse period) and the 
angular position of the pulsar are the most difficult 
parameters to determine. The period increases slowly 
with time and, in general, the pulsar moves in the sky 
(proper motion). These two parameters are discussed 
together l>ecause they are measured simultaneously. Their 
effect on the times of arrival of the pulses at a chosen 
observing point can be easily noted if the signal is 
sufficiently strong. An error in position of 1 arc-s will 
produce an arrival-time error on earth of 2.5 ms, 
ordiniTily an easily detected error. Also, apparently small 
errors in period can multiply over a large time extent so 
that predicted arrival times are seriously in error. For 
example, consider a pulsar in which the period F — 1 s 
and the uncertainty in the period is given by AP/P - 1 X 
10 ®. A prediction of an arrival time 1 week later may well 
Ik? required in this navigation scheme. In this case, the 
prediction will differ from reality by 600 ms. 

Since period and position do affect the results of time- 
of-arrival measurements, they are solved for simulta- 
neously, Note that pulsar position affects earth-bound 
measurements because the earth is in motion. The analysis 
of the results of a series of TOA arrival measurements 
requires an ephemeris of the earth in a coordinate system 
whose origin is the center of mass (barycenter) of the solar 
system. Therefore, although pulsar positions can be 
determined to an accuracy of 0.01 arc-s, these errors apply 
only within the reference frame defined by the ephemeris. 
Several observing programs have been or are being 
carried out to determine period, position, and slow 
changes in time of both quantities. The most accurate 
observing program known to this author is carried out at 
the Jet Propulsion Lalxiratory. This observing program is 
dcscril)ed theoretically in Refs. 4 and 5, and some results 
are presented in Ref. 6. 

B. Radio Spectra 

The peak flux density Bp (measured in W/m^/Hz) of a 
pulse changes with observing frequency. One generally 
sce:» a decreasing pulse strength with increasing frequency. 
A common measure of the spectrum of a pulsar is the 
energy contained within a pul.se 

m - Bp{F)if{t)di 

This ({uantity varies from pulse to pulse (see section H), 
hut a particular value is obtained after averaging over 
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AVERAGE PULSE ENERGY, 



Fig. 3 Average pulse energy versus observing frequency: (a) 
Pulsar PSR 1133 + 16. (b) PSR 0329 ^ 54 


many puLes. Of course, this average energy changes with 
frequency, and two examples of £(F) are presented in Fig. 
3. The t^havior in Fig. 3a is typical, but exceptions do 
occur (Fig. 3b). 

C. Pulse Shapes 

Imagine a simple receiver in which an RF signal is 
square-law detected, then low-pass filtered to minimize 
noise fluctuations. The output of this receiver is propor- 
tional to the average RF energy at a given time, averaged 
over a time determined by the low-pass filter. Outputs of 
such a receiver are shown in Fig. 4 for several pulsars. 
This particular sequence of pulse shapes was chosen to 
illustrate the variety of .shapes among pulsars. One should 
also note the differences in pulse width, or sharpness, and 
the different peak S/N ratios. 

D. Polarization 

The emissions of pulsars are generally randomly 
polarized, obtaining some linear polarization in .some 
cases. A good discussion of this property of several pulsars 
is given in Ref. 7. If one were forced to use a receiving 
antenna sensitive to one polarization only, it should l>e a 
circular polarization. The polarization may l>e linear or 
circular if both orthogonal monies are received. 

E. Background Radiation 

Outside the earth's atmosphere, synchrotron radiation 
from the galaxy contributes the most radiation to the 
general (and unwanted) system noise. This < ontribution Tg 
to the equivalent noise temperature Ts of the receiver 
system varies with position in the sky. At 2(K) MHz, the 
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galactic center is at a temperature of 12(K) K, while the 
galactic poles are near 13() K. Antennas used in this 
navigation system will l>e of low enough gain that the 
average temperature Tg will sufiice. The average temprra- 
turr varies with frequency, dropping from 1350 K at 100 
MHz to 25 K at 500 MHz. Using data presented in Ref. 8. 
the system temperature can l>e represented as 

Ts - T„ + Tg 
or 
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where is the ecjuivalent noise temperature of the first 
receiver amplifier and F is the observing frequency in 
megahertz. 


F. Dispersion 

The index of refraction in an ionized medium changes 
with fre(|uency. The interstellar medium consists predomi- 
nantly of ionized hydrogen with densities low enough that 
collisions arc not important. In this case, the group 
velocity of the pulse of energy is 


if f 2 ^ c is the iree-space velocity, and 8.1 X 

10*^ where F\ is in hertz and the electron density, 
is in electrons per cubic centimeter. In traveling a 
distance L (measured in parsecs, where 1 parsec = 3.08 
X 10*^ km), a time delay of AT s relative to the free- 
spac*? propagation time is introduced: 


AT 


4.1 X 10'5 



This time eWay changes with frecpiency, and fo. small 
changes in frecjuency, we find that 



1.21 X 10 



MHz/s 


(13) 


where F is now in megahertz; is the number of 

electrons in a column 1 cm^ in cross-section and /. parsec- 
long and is referred to as the dispersion measure. 


\ dispersed signal will then sweep through a receiver 
pa.ssband, as illust*-atcd in Fig. 5 The instantaneous 
bandwidth of the dispersed pulse w*!l then be less than for 
the original undispersed pulse, and distortion is the result. 
Increased dispersion produces not only increased delay 
and distortion but also lower peak S/N ratios Therefore, a 
wider bandwidth does not necessarily mean the increa.sed 
S/N ratio u.sually obtained in radiometry. Dispersion is the 
phenomenon which mosi seriously affects the signal; this 
phenomenon is dealt with directly in Section IV. 


G. Scintillation 

The shape and amplitude of individual pulses of 
radiation differ significantly from pulse to pulse for a 
particular pulsar (se<‘ e.g., Ref. 9). However, when one 
averages over several hundred pulses, the.se variations in 



TIME 



TIME 


Fig. 5 Distortion duo to disporsion in the interstellar medium 


shape and amplitude approach stable shapes and values. 
Inhomogeneities in the interstellar plasma, acting as large 
moving lenses, impose variations of several minutes on the 
pulse amplitude (see, e.g.. Ref. 10). These variations wash 
out after several hours of a' craging, and are not 
considered a deleterious effect in this study. 


IV. The Evaluation of Two Schemes 

A. Introduction 

The present objective is to choose a ni three pulsars 
which yields the lowest navigation errors. The pulsar 
period and spectrum must be known. Then, for a 
particular system temperature Tj, a peak antenna 
temperature of the pulsar, and a post-detection 
integration time Af, one can .superimpose N pulsps and 
obtain the first term of Eq. (5). Thronghoiit this analysis, it 
was assumed that ihe aiitenna gain G *= 1 and that the 
antenna is sensitive to only one polarization. The receiver 
is a direct radiometer, as descril>ed in Fig. 2, with the 
coi relation taking place on earth. With these assumptions, 
Tp can 1)0 expressed at a particular frequency as 


T,iF) 


Bp{F)A^ 

2k 


whe-e k is Boltzmann’s constant, is the effective area of 
the antenna, and Bp{F) is the o^'ak flux density of the 
pulse S’ cf r. 4rrF^Af/('^, 
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rhf wstom *'Mnprratur(' l\ i. calculated Tf^ = 300 K 

ill K(j il It is assuinctl that l\ > T^.. Svihstitution of K(j. 

14 into E(j. 5 . taking into account the superposition of 
\ pulses. \ ields 

SWh \B, ,(-(;/ XiA/^/A/')-A/' ' ^ 

where the intei»ral in K<j. .5' has In'en written in discrete 
form since the tempi ttes are available only in discrete 
form. The Munplin^ intcrsal is Af'. 

At first ^laii' c. it would seem desirable to make AF and 
AF/ as lar< 4 e a.s possible When the time il takes the 
dispersed siirnal to sweip through the passband l)ecomes 
as lar^e as the pulse width, (he sharp slniclure in the 
pulse is lost. I'nless dispersion effects are taken out. an 
unlimited HF bandwidth will result - noorer time 
resolution. Furthermore, close scnitiny of \Li\. shows 
that the expression is, except for the effects of smoothing, 
independent of AF/. Apparently, the noise contribution 
appears not so much as a S/N ratio but as the amount of 
noise eneri»\ in the time inlerxal needed for a sinnifi< anl 
decorrel ’ of the noise. This time is Af in F(j. (15>. The 
liandwidth AF/ d<H\s ha\e an effect, though, on the 
sharpness of the pulse. Limits on AF and AF/ must then 
\w found. 

In the next two sections, navigation errors are calcu- 
lated foi different comliinations of pulsars. In the Section 
B, it is assumed that the receiver cannot correct for 
dispersion in tht RF passbaiui, and this results in narrow 
baiulwidths and I. ri;er position errors. In the Section C, a 
receiver is consi meted to eompensale for dispersion 
effects, so smaller errors arc obtained by using larger 
baiulwidths. Frrors peculiar to each ptiisar were calcu- 
lated using F(j. t5.'. Position errors obtained using a given 
set of three pulsars were determined as de.scrilied in 
Section nC. 

^ B. The Dispersionless Receiver 

The simplest spacecraft receiver possible is considered 
first to provide insight into the navigation problem. The 
receiver is represented in Fig. 6. The antenna eonfigura- 
tion remains unspecified, but il is assumed to cover the 
s c elestial sj liere, each pulsar experiencing an antenna gain 

of unity. Three different receiver channels are needed 
since the spectra an 1 dispersion measures differ amonc 
pulsars. Integration or superposition of the pulse.s o\er <\ 
hiiig lime inteiA’al (hours) will i>e necessary, so the output 



of cMch receiver channel is fed into its own digital 
sai ling and storage device. Sampling and the storage 
addr^'s^ are coordinated and controlled b a controlled 
oscillator. Since the spacecraft is in motion, significant 
doppler shifts will lie experienced. Failure to adjust the 
sampling rale so that the same numl>er of samples are 
taken during each pulse period will re.siilt in severe 
smearing of the pulse. Therefore doppler shift infor.^iation 
IS sent to the spacecraft on the uplink l>efore an observing 
session is to l>eg»n. The resulting integration is sent to 
earth-based computers vl • the downlink along wiln the 
lime at which .sampling started. The numl>er of pulses 
superimposed depends on how rapidly the expected 
doppler .shift is changing. One could, in principle, 
superimpose pulses for several days if no changes were 
occurring. During encounters with planets, superposition- 
ing may last only I h lx?fore an “update" of the controlled 
o.sciliator is necessary. Further superposition then takes 
place on earth. A monitor of the fijquenc*y standard 
driving the controlled oscillator is maintained. Slow drifts 
in the standard are ccounted for in the information 
supplied to the controlled oscillator. 

The performance of this system was evaluated in two 
st.,ges. First, the error ^ was calculated at several 
receiver fre<piencies I)etwe< n 100 and 1650 MHz for the 
first 18 enlrie.s in Table 1. The last nine entries are not 
ol >cr\able at Goldstone. The .spectra of these pulsars are 
know'll through a compilation of data taken from the 
literature on pulsars. The compilation presented here does 
not appear in the literature. The spectra were approxi- 
mated by power-law functions within oi^e or more 
segments of the 80- to 10,000-MHz frequency range. 
These segments are denoted in Table 1, the last segment 
in each row ending at 10,000 MHz, The low-frequency 
limit F/^ and the appropriate exponent a are listed for 
each .segment, in segment 1, the average pulse energy is 
given l)v /r(F///F)“, where F// is the highest frequency in 
segment 1 (300 MHz in the case of pulsar 0628, e.g.), and 
7, is the energy at F//. In segments 2 and 3, the energy is 
given by /,(F/ /F)®. Although more than eigh^v pulsars are 
known, only the strongest are considered here. 

System noise temperatures were calculated from Eq. 
(12). The integrals were calculated from templates 

measured at 2388 MHz in a JPL research program on the 
first 18 pulsars in Table 1. Although the templates y(0 are 
slow functions of frequency, it was assumed that they do 
not vary significantly with frequency l>ecause of a lack of 
information at lower frequencies. The RF bandwidth AF 
was slowly widened until f{t) was significantly smewthed. 
The appropriate smoothing time was calculated using the 
dispersion mea.su res peculiar to a given pulsar. Also, the 
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Fig. 6 The disporsionless receiver 


post-detection bandwidth was slowly decreased until fit) 
was affected. In this way, limits on AF and AF/. were 
placed on each receiver ( hannel for several possible 
operatiNg frecjuencies. It was assumed that pulses were 
superimposed for 21 h. This completed the first stage. 

Trie pulsars were lirought together three at a time in 
the second stage of the analysis. The positions u.scd are 
readily found in the literature. Eigenvalues of the inverse 
covariance matrix were calculated for \hc case F^ = 
4, corresponding to solving >r the 95%-coniidencc 
ellipsoicl (see the discussion of E<p 9), A total of 816 
combinations were anal v zed. The results were considered 
in order of increasing ellipsoid volume. The lowest errors 
obtainable w ith this system are presented in Table 2 along 
with the three pulsars used to obtain the spacecraft 
position, lamiting values of A/- and AF; are presented 
with expected arrival-time errors and S/N ratios. The size 
of each princ ipal axis is compared to the diameter Dj of 
Jupiter (1.5 X lO'*' km). 

'Fhe inade(|uacy of this scheme is readily apparent in 
Table 2, since the conditions for the validity of Ecj. (5) are 
not met. The values of gt^ are greater that; the pulse width 
in each case. This is also exhibited by S/N ratios less than 
1.0. \ more promising scheme is analyzed in Section C. 


C. The DisperslonXorrecting Receiver 

A large limitation in the previous scheme is the 
confinement to narrow bandwidths. It was therefore 
nostulated ihat a dispersive filter can l>e built for a 
spacecraft which will compensate for the interstellar 
dispersion. Such a system is represented in Fig. 7. The 
dispersive filter is conceptually represented as a parallel 
set of filters of width AF, followed by delay Tj. (Note that 
AT’ — 2AF,.) The progre.ssion of delays depends on the 
dispersion mea.su re of the particular pulsar. The rest of 
the receiver is identical to that previously described. 

The last nine pulsars in Table 1 were included in the 
analysis. It was hoped that inclusion of the.se southernmost 
pulsars would produce an error ellipsoid more spherical in 
shape than that presented in Table 2. Little is known 
about the .spectrum and pulse shape of these particular 
objects. The spectra were therefore conservatively 
modeled about published average strengths at 400 MHz. 
Below 4(X) MHz, the spectra were assumed to l>e flat (a * 
0.0). which means that the peak flux density at 100 MHz 
may have l>een underestimated by as much as 10. The 
detailed pulse shapes are practically unknown, so 
templates for pulsars of approximately the same period 
and overall pulse width were assigned to these less-studied 
pulsars. 
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Tabit 1. Modalt of th# tpoctra of tha ttrongatf pultart 


Entry 



1 



2 


3 

J / mVH / 

^ MHz 

a 

K , . .MHz 

a 

/•,.MHz 

a 

1 

0031 

7.0 X U ) -- 

80 

1.9 





2 

0329 

1.0 X 10-*^ 

100 

-1.3 

300 

0.0 

800 

2.75 

3 

0525 

1.0 X 10- 

100 

1.8 





4 

06lS 

1.0 X 10 --C 

80 

0.0 

300 

1.9 



5 

0823 

1.0 X 10-27 

80 

0.0 

500 

1.85 



6 

0950 

4.0 X 10-2S 

80 

1,0 

1000 

2.5 



7 

1133 

6.0 X 10 

100 

0.8 

1400 

2.1 



8 

1237 

1.0 X 10-27 

100 

0.9 

300 

1.3 



9 

1642 

1.0 X 10-2« 

80 

0.0 

150 

2.2 


1 

10 

1706 

5.0 X lO--** 

80 

0.0 

1000 

2.0 



11 

1749 

6.0 X 10 - 2 ‘ 

80 

00 

500 

2.85 



12 

1818 

1.0 X 10-27 

80 

O.o 

400 

2.6 



13 

1911 

2.0 X 10-27 

80 

0.0 

300 

2.5 



14 

1929 

4.0 X 10-2^ 

100 

0.0 

1000 

1.9 


1 

1 

j 

15 

1933 

4,0 X 10-2'' 

100 

0.0 

1000 

3.1 


t 

I 

16 

2021 

1.0 X 10-27 

100 

0.0 

1000 

1.6 


i 

j 

17 

2045 

8.0 X 10-2'» 

100 

0.0 

1000 

2.6 


1 

18 

2218 

1.0 X 10-2C 

100 

2.2 




t 

19 

0254 

5.0 X 10-27 

100 

0.0 

400 

2,0 


I 

20 

0940 

1.0 X 10-27 

100 

0.0 

40C 

2.0 


i 

j 

21 

0959 

1.0 X 10-27 

100 

0.0 

40'. 

2.0 


1 

22 

1359 

5.0 X lo ^’’ 

100 

0.0 

40( 

2,0 


i 

23 

1426 

1.0 X 10-27 

100 

0.0 

400 

2.0 


i 

i 

24 

1451 

2.0 X 10-27 

100 

0.0 

00 

1.8 


( 

1 


25 

1530 

1.0 X 10-27 

100 

0.0 

400 

2.0 

26 

1727 

1.5 X 10-27 

100 

1.0 

400 

2.0 

27 

1747 

1.0 X 10 27 

100 

1.7 

400 

2.0 


The bandwidth AF was fixed at 200 MHz and centered 
about Fq » 200 MHz. Arrival-time errors were calculated 
for 27 pulsars. A range of post*detection handwidths were 
considered for the purpose of finding limits on AF^. No 
dispersion smearing was considered since the dispersion- 
correcting filter elimio'’ktes this effect. The pulsars were 
then analyzed in > of three (2925 combinations) to 
determine navigation errors. The resulis were ranked 


according to increasing volume of the ellipsoid. The three 
best combinations are presented in Table 3. 


The analysis was repeated with a center frequency of 
250 MHz to determine the sensitivity of the errors to 
operating frequency. The results, also presented in Table 
3, show that a center frequency of 200 MHz is better than 
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T»bl« 2. Navigation arrort *or tht ditparslonlats racaivtr at 
ISO MHz. with 24 h of inttgration 


Pulsar 

•IF. 

Hz 

iF,,. 

Uz 

ms 

S/N 

ratio 

0950 

9.0 X 10' 

300 

41 

0.11 

1642 

5.6 X 10* 

380 

25 

0.12 

1749 

5.7 X 10* 

280 

66 

0.06 


250 MHz am] imply that a slightly narrower bandwidth 
may be better. 

The results presented in Table 3 indicate that with this 
low-cost navigation scheme, one can navigate among the 
outer planets with an accuracy comparable to the pre.^^ent 
navigation scheme. The system has the following charac- 
teristics: 

(1) Each pulsar is viewed with an antenna gain of 1.0. 

(2) The antenna is sensitive to one polarization. 

(3) The receiver is a direct radiometer. 

(4) Dispersion effects in the receiver passband are 
eliminated. 

(5) The average pulse shapes are known, the detection 
process is square-law, and the time-of-arrival is 
determined by correlating the detected signal with 
the average template. 


Several alterations that will improve this system are 
discussed in Section V. 


V. Concluding Remarks 

The expected navigation errors listed in Table 3 would 
be extremely attractive if they were smaller by a factor of 
10. This decrease can be most easily obtained by 
increasing the antenna gain. A broadband or frequency- 
independent antenna is needed to cover the frequency 
band of 100 to 300 MHz. Imagine, then, a trapezoidal wire 
antenna supported by two booms. Another trapezoidal 
wire antenna is rotated 90 deg in space and physically 
combined with the other wire antenna to form a 
truncated-pyramid shape. The base would be 1.5 to 2.0 m 
on a side and the “top** or small end of the pyramid about 
0.5 m. The height of the tnmcated pyramid would be 
about 0.5 m. (A 45-deg angle exists between the centeiline 
and each side of the antenna; see Ref. 11.) 

An antemia such as this would have a gain of 7 to 8 dB 
and a circular beam of 66-deg half-power width. Nine 
such antennas would be needed to cover the entire sky. 
This is probably not possible since such a structure of 
wires would interfere with the communications antennas. 
Also, complete sky coverage allows the sun's radiation to 
cause an increase in the equivalent noise temperature of 
the receiving system. Only three antennas are needed if 



Fig. 7 Tht ditptrtion-corrtcting rtctivtr 
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Table 3. Navigation errors for the dispersion-correcting receiver at Fq « 200 arj 250 MHa. AF « 

200 MHz, with 24 h of integration 


Combination 

Pulsar 


<7y, ms 

PriiK-ipal 

axis 

Si/t*. kill 

11/ 

K„ 200 

2.50 

200 

l-„ 250 


UU2 

500 

0.40 

0.74 

1 

650 

1160 

1 

1749 

600 

0 90 

1.07 

o 

IStX) 

4140 


221S 

■1«0() 

O.SO 

2.11 

3 

1520 

2130 


1()42 

500 

0.40 

0 74 

1 

660 

1230 

0 

1911 

30(K) 

l.SO 

2.23 

2 

3800 

5800 


22 IS 

4600 

0.80 

2.11 

3 

1.540 

3190 


0950 

5000 

1.80 

3.03 

1 

660 

1220 


1642 

500 

0.40 

0.74 

2 

3820 

6980 


221S 

4600 

0.80 

2.11 

3 

1590 

3900 


the spacecraft maintains a constant orientation relative to 
the star. 

This antenna system would be sensitive to both 
polarizations, so an improvement of 10 dB can be 
expected. Using combination 1 of Table 3, the low errors 
of Table 4 can l>e obtained. 

A few remarks concerning ephemeris and timing errors 
are necessary. There is always some uncertainty in the 
position of a particular pulsar This causes some uncer- 
tainty in referring the arrival time of a particular pulse at 
the earth to the origin of ♦he coordinate system. If the 
position error is 0.01 arc-s, the maximum error in arrival- 
time is 24 ^s, a barely noticeable effect. The spacecraft 
sees the same position error. At the distance of Pluto, the 
maximum error in referring an arrival time to the origin 
(or alternatively, in the solutions to Eq. 4) is magnified to 
alxuit 0.9 ms. This is a significant systematic error. 

Ephemeris errors have an effect on this scheme in that 
’.he position of earth relative to the origin of the 


Table 4. Navigation errors for F(> ■ 200 MHz, AF * 200 
MHz, and G « 10, with 24 h of integration 


Principal axis 

Size, km 

1 

65 

2 

180 

3 

150 


coordinate system (usually the center-of-mass of the solar 
system) is in error by about 35 km in current ephemerides. 
This amounts to an error in arrival time of about 120 /is. 
Although this is a small error for the navigation system, it 
is not small for earth-based arrival time measurements; 
i.e., pulsar arrival-time data could be used to improve an 
ephemeris of earth. 

The determination of the pulsar models is not very 
sensitive to the location of the receiving antennas. A 1-km 
error in station location causes a 3.3 /ts error in arrival 
time. This error is not significant in current measurements. 

Clock errors affect arrival time measurements directly. 
These errors should be 1/10 or less of the desired arrival 
time accuracy. This applies to spacecraft and earth-based 
clocks. 

The definition of epoch appears in this navigation 
scheme as well as in the schemes currently in use. During 
24 h of integration, a spacecraft can move as much as 10^ 
km. One experiences the same conceptual problem in 
current systems when integrating range and doppler shift 
errors over 12 h of spacecraft tracking. 

The acceleration of the spacecraft creates a need for 
second-order doppler shift corrections. If the spacecraft is 
accelerating uniformly, the actual phase of the pulse train 
drifts relative to the phase of the zero-order pulse train 
(constant period modified by a constant doppler shift) by 
an amount proportional to the square of time. Represent- 
ing the velocity in the direction of the pulsar by uo + 
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I 


I 

i 

- to), the phase drift, expressed as a fractional part of the have been observed (see, e,g.. Ref. 12). Although no such 1 

rest-frame period, becomes large jumps have been observed in the pulsars included in i 

this study, such possibilities caniiol be eliminated. In 
~~ /} -f. £2.\ ^ addition, pulsar models of period and position need to be 

2c Pq \ c / updated because of small errors in the models caused by 

errors in arrival time measurements. 

where c is the velocity of light and Vq is the zero-order 

velocity. The integrated pulse is smeared over an entire Finally, it should be pointed out that there is an 
pulse period if (1) the polynomial of Fig. 6 is zero-order. ambiguity in this navigation technique. The measurement i 

(2) the integration time t - to is 24 h, and the acceleration of the pulse arrival time for a given pulsar relative to a 

is only 0.08 m/s^. Therefore, acceleration must be given epoch is the measurement of the phase of a pulse 

allowed for in the polynomial generator. It may be train. There is an ambiguity of an integral number of pulse 

unreasonable to expect a uniform acceleration, so a periods NP in the arrival time. In fact, we determine a 

second-order model of the velocity may be needed. series of planes separated in space by Pr km (3 X 10^ km 

for a 1-s pulsar). In using three pulsars, we then determine 
This navigation scheme does require a continuation of a three-dimensional grid of possible spacecraft locations, 

earth-based measurements of pulsar arrival times. Discon- Simultaneous use of current navigation techniques would 

tinuities or sudden jumps in the period of two fast pulsars eliminate these ambiguities. 
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